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Abstract. We present an intrinsic formulation of the kinematic prob- 
lem of two n— dimensional manifolds rolling one on another without twist- 
ing or slipping. We determine the configuration space of the system, 
which is an "( T ^+ 3 ) —dimensional manifold. The conditions of no- twisting 
and no-slipping are decoded by means of a distribution of rank n. We 
compare the intrinsic point of view versus the extrinsic one. We also show 
that the kinematic system of rolling the n-dimensional sphere over R™ 
is controllable. In contrast with this, we show that in the case of SE(3) 
rolling over se(3) the system is not controllable, since the configuration 
space of dimension 27 is foliated by submanifolds of dimension 12. 



1. Introduction 

Rolling of surfaces without slipping or twisting is one of the classical 
kinematic problems that in recent years has again attracted the attention of 
mathematicians due to its geometric and analytic richness. The kinematic 
conditions of rolling without slipping or twisting are described by means of 
motion on a configuration space being tangential to a smooth sub-bundle 
that we call a distribution. The precise definition of the mentioned motion in 
the case of two n-dimensional manifolds imbedded in R m , given for example 
in [11], involves studying the behavior of the tangent bundles of the manifolds 
and the normal bundles induced by the imbeddings. This approach leads to 
significant simplifications, for instance, it suffices to study the case in which 
the still manifold is the n— dimensional Euclidean space. The drawback is 
that the geometric descriptions depend strongly on the imbedding under 
consideration. 
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However, so far little attempts have been made to formulate this prob- 
lem intrinsically. An early enlightening formulation is given in [2], in which 
the authors study the case of two abstract surfaces rolling in the above de- 
scribed manner. This is achieved by means of an intrinsic version of the 
moving frame method of Elie Cartan which, for this case, coincides with 
the classical intrinsic study of surfaces, see [12]. One of important results 
established in [2] is the non-integrability property of a rank two distribution 
corresponding to no-twisting and no-slipping restrictions, namely, if the two 
surfaces have different Gaussian curvature, then the distribution is com- 
pletely non-integrable and moreover is of Cartan-type, see [I]. A control 
theoretic approach to the same problem, studied in [TJ, has the advantage 
that the kinematic restrictions are written explicitly as vector fields on ap- 
propriate bundles. 

We present a generalization of the kinematic problem for two n— dimensional 
abstract manifolds rolling without twisting or slipping via an intrinsic for- 
mulation. We define the configuration space of the system, which is an 
^Y^— dimensional manifold and which is a direct analogue to the one 
found in the references [2] and [I]. We give several equivalent definitions 
of rolling motion involving intrinsic characteristics and those that depend 
only on imbedding and discuss their relations. This new definitions permit 
to determine the imbedding-independent information contained in the ex- 
trinsic definition of the rolling bodies problem presented in [H] . Moreover, 
we relax the smoothness condition of the rolling map up to absolutely con- 
tinuity. This allows to enlarge the class of mappings under consideration, 
still giving the possibility to apply the fundamental theorems of differential 
geometry and control theory without changing drastically the main classical 
ideas of rolling maps. The conditions of no-twisting and no-slipping define 
a distribution of rank n in the tangent bundle of the configuration space. 
We write explicitly the distribution as a local span of vector fields defined 
on the configuration space. We test the bracket generating condition of the 
above mentioned distribution on the known example [H] of rolling the n- 
dimensional sphere over the n-dimensional Euclidean space and the special 
group of Euclidean rigid motions SE(3) rolling over sc(3). As a result we 
obtain the controllability of the first system and the non controllability of 
the latter. 

The structure of the present paper is the following. Section [2] is an in- 
troductory section where we collect necessary definitions and discuss the 
motivation for the reformulations of kinematic conditions of no-twisting and 
no-slipping for the rolling problem. We present two formulations and show 
their equivalence. Section 3 gives a good starting point for comparing differ- 
ent approaches, known in the literature for 2-dimensional rolling manifolds. 
In Section 4 we give the main formulation of extrinsic rolling as a curve on a 
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configuration space defined as a direct sum of principal bundles over Carte- 
sian product of two rolling manifolds and we prove the equivalence of new 
extrinsic definition of rolling with the previous ones and deduce the intrinsic 
definition of a rolling map. We also prove some theorem enlightening the 
imbedding independent information contained in the principal definition of 
the extrinsic rolling. Section 5 is dedicated to the construction of two distri- 
butions in the tangent bundle of the configuration space. These distributions 
encode the no-twisting and no-slipping kinematic conditions of extrinsic and 
intrinsic rolling. The rolling, both extrinsic and intrinsic, can be written as 
a curves in configuration space tangent to the corresponding distributions. 
In Sections 6 and 7 we present detailed calculations for the two aforemen- 
tioned examples: rolling the n- dimensional sphere over the n-dimensional 
Euclidean space and rolling SE(3) over se(3). In the first case the distri- 
bution is bracket generating, coinciding with the result obtained in [JJT In 
the second case we obtain that the configuration space, of dimension 27, is 
foliated by 12 dimensional sub manifolds. 

2. Definition of rolling map for manifolds imbedded in 

Euclidean space 

2.1. Rolling without twisting or slipping for imbedded manifolds. 

We start from the classical definition of rolling without slipping or twisting 
of one manifold over another manifold inside the Euclidean space. 

Let us start with some notations. Throughout this paper, M and M will 
always be oriented connected Riemannian manifolds of dimension n. By the 
well known result of Nash, see [H], there are isometric imbeddings of M and 
M, denoted by i and T respectively, into IR n+t/ for an appropriate choice of v. 
Here and in what follows M. n+U will always be equipped with the standard 
Euclidean metric and standard orientation. As long as there is no possibility 
for confusion, we will identify the abstract manifolds M and M with their 
images under the corresponding imbeddings. The imbedding of M into W l+U 
splits the tangent space of ~R n+u into a direct sum: 

(1) T x R n+u = T X M © T X M L , x e M. 

In general, any objects (points, curves, . . . ) related to the manifold M will 
be marked by a hat (~) on top, objects related to M will be free of it, 
while terms related to the ambient ~R n+u space carry a bar ( ~ ) . We use 
Isom(M) for the group of isometries of M, and Isom + (M) for the group of 
sense preserving isometries. 

We start by given the definition of rolling without twisting and slipping 
as found in |llj . 

Definition 0. Let M,M be submanifolds of M. n+l/ . Then, a differentiable 
map g : [0, t] — > Isom(IR n+I/ ) satisfying the following conditions for any 
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t G [0, r] is called a rolling M on M without slipping or twisting. The 
rolling conditions: 

• there is a piecewise smooth curve x : [0, r] — > M, such that 

— g(t)x(t) EM, _ 

— T g{t)x{t) (g(t)M) = Tg {tMt) M. 

• Furthermore, the curve x(t) := g(t)x(t) satisfies the following 

— no-slip condition: 

g^gity^it) = 0, 

— no-twist condition, tangential part: 

d(g(t)g{ty l )T m M C T Q {g{t)g{t)- x M) L , 

— No-twist condition, normal part: 

digitMty^T^M 1 C T Q {g{t)g{t)- l M) 

Remark 1. In the previous definition, we explicitly state that g : [0,r] — > 
Isom(IR n+i/ ) is differentiable. This is not stated in [11], but conditions con- 
taining g are required to hold for all t. Also, a minor inaccuracy in the 
no-twisting conditions is corrected. 

It is clear that Definition is of extrinsic nature. Thus, in order to 
obtain an intrinsic formulation of the rolling problem, we want to change 
the original definition as follows: 

(1) Making x(t) part of the data of the rolling: The reason is to give a 
local character to conditions of rolling without twisting or slipping. This 
will emphasize the dependence of the rolling not just on the isometry g but 
also on a curve x along which the rolling of M on M can be realized. In 
some particular cases, this may lead to small changes in terminology. The 
following example illustrates these ideas. 

Example 1. Consider the submanifolds of M 3 , defined by 

M= {(x!, 1- cos 9, sinO) Gl 3 | X\ G E, 6e [0,2vr)}, 

M = {(xi,x 2 ,0) G M 3 | Xi,x 2 GM, }. 

The rolling map 

x 2 M-> I (x 2 — 1) cost + x 3 sin t + t + cos 2t 
x 3 J \ — (x 2 — 1) sin t + x 3 cost — sin 2t 

describes the rolling of the infinite cylinder M on M along the avaxis with 
constant speed 1. Then there is an infinite choice of curves x(t) G M, given 
by 

x(t) — (xi, 1 — cost, sin t), Xi G R 




ROLLING MANIFOLDS 



5 



along which the rolling g can be realized. However, if we make x(t) as part 
of the data, then each choice of the curve x(t) will correspond to different 
rollings (x(t),g(t)) . 

(2) Relaxing the differentiability conditions for g(t): We think that the con- 
ditions of differentiability of g(t) for all t G [0,r] and piecewise smooth- 
ness of x(t) are too restrictive. The requirement that (x, g) : [0, t] — > 
M x Isom(M n+!/ ) is absolutely continuous or Lipschitz seems more natu- 
ral, since this allows us to implement results from control theory, see Sub- 
section [3T2l In this context, absolute continuity of a curve on 
M x Isom(M n+ly ) is considered with respect to the parameter t, as in [TJ 
Chapter 2]. 

(3) Introducing orientability assumptions: In order to have a connected con- 
figuration space, we exploit the orientability assumption of M and M. Since, 
as mentioned before, the rolling conditions will be local, we may choose an 
orientable neighborhood of the starting point even on any non-orientable 
manifold. We will use this to impose some practical restrictions to the defi- 
nition of a rolling. 

• Since g(t) is continuous, it is either always orientation preserving or 
orientation reversing isometry of IR n+t ' for all t. Given a rolling g(t) of 
M on M, we may assume that g(t) is always orientation preserving by 
changing the orientation of M n+y . To obtain an orientation preserving 
rolling from an orientation reversing rolling g(t) of M on M, pick any 
constant orientation reversing isometry go of M n+i/ . Then gog{t) is 
an orientation preserving rolling of M on go(M). 

• It is intuitively clear that for a fixed t, d x ^g(t) maps elements from 

T x (t)M to T x [t)M and elements from T^M 1 - to T^M 1 - (for more 
details see Subsection 12.21) . Hence, the matrix form of d x (t)g{t) splits 
in the following way: 



T x[t) M T x{t) M L 

A(t) \ T x{t) M 



wyv ; V ) T x{t) M A 

Since g(t) is orientation preserving, both d x ^g(t)\T x{t) M and d x (t)g{t)\ T ^ 
are either orientation preserving or orientation reversing. By conti- 
nuity, d x (t)g{t)\T x(t) M is either orientation preserving or orientation 
reversing for all t. We will require that d x ( t )g(t)\T x(t) M is always ori- 
entation preserving. If d x ( t )g(t)\T x{t) M is orientation reversing, pick 
any constant orientation preserving isometry go : M. n+U — > W n+L/ so 
that 

dg \ TU : TM -»■ T(g M) 



.r(f) 



A I- 
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is orientation reversing. It is sufficient to show that it reverses the 
orientation at one point in order to show that it reverses orientation 
at all points due to the fact that M is oriented. Then gog(t) will be 
a rolling of M on g M which is orientation preserving on the tangent 
space at x(t). 

Implementing the above changes to Definition 0, we obtain the following, 
from which several equivalent reformulations will be presented later. 

Definition 1. A rolling of M on M without twisting or slipping is an ab- 
solutely continuous curve (x,g) : [0,r] — > M x Isom (l n+!/ ), satisfying the 
following conditions: 

(i) x{t) := g{t)x{t) G M, for all t e [0,r]. 

(ii) T m (g(t)M) = T m M, for all t e [0,r]. 

(iii) No slip condition: g(t) o g~ l (t)x{t) = 0, for almost every t. 

(iv) No twist condition (tangential part): 

d(g(t) o g~\t))(T x(t) M) C T (g(t) o g^^M) 1 , 

for almost every t. 

(v) No twist condition (normal part): 

d(g(t) o g-\t))(T m M^) C T (g(t) o g-\t)M), 
for almost every t. 

(vi) d x{t )g(t)\ Tx(t) M ■ T x[t) M ->■ T x{t) M is orientation preserving, for all 
tE[0,r]. 

We omit, from now on, the words "without twisting or slipping", just 
writing "a rolling of M on M". Furthermore, for given curves x(t) and x(t) 
in M and M, respectively, the expression "a rolling of M on M along x(t) 
and x(t) v will mean a rolling (x,g) : [0, r] — > M x Isom + (lR n+ ^) so that 
g{t)x{t) = x{t). 

Remark 2. The definitions we will be working with ignore physical restric- 
tions given by the actual shapes of the manifolds. Intuitively, if we think 
of the manifolds in Definition [T] as physically touching along the curves x(t) 
and x(t) and rolling according to the isometry g(t), then we cannot rule out 
the possibility that there might be non-tangential intersections between the 
manifolds other than the contact points. 

Example 2. Consider the imbedded surface 

M = {(x 1 ,x 2 ,x 3 ) e M 3 | x\ - x\ + x 3 = 0, x\ + x\ < 1}, 

and M = M 2 , imbedded as an affine plane. Assume that both manifolds M 
and M carried the induced metric. We can clearly define a rolling of M on 
M in terms of Definition [JJ but there is no way to connect the saddle point 
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in M with any point in M without there being intersections between the 
surfaces. 

2.2. First reformulation. We aim to give a definition of the rolling of 
M on M in a way that is more fruitful for future considerations. We fix 
some notations first. According to the splitting (JTJ, any vector v G T x M. n+ ", 
x G M, can be written uniquely as the sum v = v T + f -1 , where v T G T X M 
is tangent to M at x, while v 1 - G T X M^ is normal. Analogous projections 
can be defined for M. 

Let V denote the Levi-Civita connection on M or on M. The context 
will indicate on which manifold the connection is defined. The "ambient" 
Levi-Civita connection on ~R n+u is denoted by V. Note that if X and Y are 
tangent vector fields on M, then 

V x Y(x) = (VxY(x)) T , xeM, 

where X and Y are any local extensions to ~R n+u of the vector fields X and 
Y, respectively. Similarly, if T is a normal vector field on M and X is a 
tangent vector field on M, then the normal connection is defined by 

ViT(x) = (VxT(x)) ± , xeM, 

where T is any local extension to ~R n+u of the vector field T. Equivalent 
statements hold for M. If no confusions arise, we will use capital Latin 
letters X, Y, Z to denote tangent vector fields and capital Greek letters T, 
for notation of normal vector fields. 

For a fixed value of x G M and a fixed vector field Y, the vector VxY(x) 
only depends on the value of X{x). Therefore, for v G T X M, we will use V V Y 
or V v Y(x) to mean VxY(x), where X is an arbitrary vector field satisfying 
X(x) = v. We will use the same convention when V is interchanged with 
V-K 

If Z{t) is a vector field along x(t), we will use j[Z(t) to denote the covariant 
derivative (corresponding to V) of Z(t) along x(t), and for any normal vector 
field ty(t) along x(t), ^-^(t) denotes the normal covariant derivative (see [TJ 
p. 119]). Recall that if M is imbedded isometrically into M n+I/ , then 

where Z(t) and ^(t) are tangential and normal vector fields, respectively, 
along a curve in M. 

We say that a tangent vector Y(t) along an absolutely continuous curve 
x(t) is parallel if j^Z(t) = for almost every t. Notice that it is possible to 
define the notion of parallel transport even though the derivative x(t) exists 
only almost everywhere, see, e. g., Existence and Uniqueness Theorem in [Till 
Appendix C]. Namely, for any absolutely continuous curve x : [0,r] — > M 
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and for any v G T x ( to )M, < t < r, there exists a unique absolutely 
continuous tangent vector field Z(t) along x(t), such that Z(t) is parallel 
and satisfies Z(t ) = v. 

We say that a normal vector field ^{t) along x(t) is normal parallel if 
^-(t)\& = for almost every t. A normal analogue of parallel transport is 
defined likewise. 

We are now ready to give a new formulation of the rolling map. 

Definition 2. A rolling of M on M without slipping or twisting is an ab- 
solutely continuous curve (x,g) : [0, r] — > M x Isom + (M n+v ) satisfying the 
following conditions: 

(i') x(t) := g(t)x(t) G M, 

(ii') dg(t)T x{t) M = T m M, 

(iii') No slip condition: x(t) = dg(t)x(t), for almost every t. 

(iv') No twist condition (tangential part): 

dg(t)j t Z(t) = jdg{t)Z, 

for any tangent vector field Z(t) along x(t) and almost every t. 
(v') No twist condition (normal part) : 

for any normal vector field *&(t) along x(t) and almost every t. 
(vi') d x (t)g{t)\ Tx{t)M ■ T x (t)M ->■ T^M is orientation preserving. 

Lemma 1. Definitions^ and\^ are equivalent. 

Proof. Since (i) and (i') coincide, we begin by proving the equivalence of (ii) 
and (ii'). Restricting the action of g(t) to M, we observe that the differential 
d x (t)9(t) maps T x mM into T 9 m x m (g(t)M) by definition, and hence (ii) holds 
if and only if (ii') holds. 

In order to prove the equivalence between (iii) and (iii') we write a curve 
g(t) in Isom + (M n+l ') as 

g(t):x\-> A(t)x + f(t), xeR n+ ", 

where A : [0, r] -> SO(n + v) and f : [0, r] ->■ R n+V . Thus d x g(t)v = A(t)v, 
v G T x R n+u , and we get 

git) o g-\t) x(t) = g(t)x(t) = A(t)x(t) + f{t) 

= j t (A(i)x{t) + f{t)) - A{t)x{t) = - dg{t)x{t). 

whenever x(t) is defined. Hence g(t) o g~ l {t)x{t) = almost everywhere if 
and only if x(t) = dg(t)x(t) almost everywhere. 
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Before we continue with the final two conditions, notice that (ii') actually 
states that both dg{t){T x[t) M) = T m M and dg{t){T x{t) M^) = T^M 1 - hold 
due to the splitting ([1]). Hence, the inverse differential dg~ l {t) = (dgit))^ 1 
also maps tangent vectors to tangent vectors and normal vectors to normal 
vectors. This allows us to restate (iv) and (v) as the conditions 

(dg(t)v T ) T = 0, and (dg(t)v ± ) ± = 0, 

holding for almost every t and for any v G T x ^M. n+I/ , decomposed as the 
sum of v T G T x (t)M and v 1 G T^M 1 - via the splitting (Op). We calculate 

= (dg(t) Z(t)) T = (± (dg(t) Z(t)) - dg(t) (^Z(t)j ) T 

= 2- t dg(t)Z(t)-dgj t Z{t) 

for any tangent vector field Z{t) along x(t), for any value of t where x(t) is 
defined. By similar calculations, using a normal vector field along x(t), 
we obtain 



□ 



Remark 3. The following observations are useful for the understanding of 
the nature of a rolling map. 

• The proof of Lemma ED shows that indeed condition (ii') is equivalent 
to the statement 

dg{t)T x{t) M ± = T m M ± . 

• Condition (iv') is equivalent to the requirement that any tangent 
vector field Z(t) is parallel along x(t) if and only if dg(t)Z(t) is 
parallel along x{t). As a consequence, this condition is automatically 
satisfied in the case of one dimensional manifolds. 

• We can reformulate (v') in terms of normal parallel vector fields. 
Namely, condition (v') is equivalent to the statement that any normal 
vector field is normal parallel along x(t) if and only if dg(t)^>(t) 
is normal parallel vector field along x(t). Thus, if the manifolds 
are imbedded into Euclidean space and the codimension is one (i.e. 
v = 1), condition (v') always holds. 



3. Previous intrinsic descriptions of rolling maps dimension 2 

The aim of this Section is to present the different intrinsic formulations 
of a rolling map appearing in literature for two dimensional manifolds. The 
two best known formulations are given in fTj [2] . We start by introducing the 
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configuration space of the rolling for the general case of n dimensional mani- 
folds and then proceed to describe the previously mentioned two dimensional 
situation. 

3.1. Frame bundles and bundles of isometries. Let M and M be ori- 
ented connected Riemannian manifolds of dimension n. We introduce the 
configuration space Q of the rolling, which can be thought of as all the 
relative positions in which M can be tangent to M. Define the principal 
SO(n)-bundle over MxMby 

Q= \qelsom+(T x M,T $ M) 



x G M,x G M 

Here Isom^fV, V) denotes the group of linear orientation preserving isome- 
tries between the oriented inner product spaces V and V. 

The principal SO(n)-bundle structure of the configuration space Q can 
be also described in the following way. Let F and F be oriented frame 
bundles of M and M, respectively, with the obvious principal SO(n)-bundle 
structures. Consider F x F as a bundle over M x M with SO(n) acting 
diagonally on the fibers. Then, we can identify Q with (^F x Fj / SO(n) by 
the following map. Let {e J (x)}" =1 be a frame in F at x G M and similarly 
let {e-i(x)}f =1 be a frame in F at x G M. To the equivalence class 

(fe(x)}^ =1 ,{e,(x)}r =1 )-SOH 

we associate the unique isometry q G lsom.Q (T X M,T$M) satisfying 
(2) ei(x) = qei(x), i = l,...n. 

Clearly, this construction does not depend on the choice of a representative 
of an equivalence class of (^F x Fj / SO(n). Conversely, given an isometry, 

there exists a unique equivalence class of frames satisfying (121). 

The left and right action on fibers of Q is induced by the inverse left 
action on F and left action on F, respectively. More precisely, an element 
Aq G SO(n) acts on an isometry q G Q from the right or left, giving q\ = qAo 
and q2 = A q, respectively. The isometries q\ and q2 are defined by 

ei(x) = qi A^dix), A ei(x) = q 2 e^x), 
where {{ej(x)}j =1 , {ej(x)}™ =1 ) is any basis satisfying (j2|). 

Furthermore, let [{fj(x)}j = i, {fi(x)}i=ij G (F x F)]^) be any other 
pair of frames and the matrix representation A = (a^)" . =1 , of an isometry 
q is given by 

A= (a^. =1 = (j*^)qfM)) n i]=1 ■= ((f^),qf)(x))Y ij=1 e SO(n), 
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where /* stands for the 1-form dual to the vector field /j. Then, in this 
basis, 

(fi(x)qifj{x)) n =AA , and (f*(x) q 2 J,(x))" = A A. 

\ / i,j=l \ / i,j=l 

Since Q is a principal SO(n)-bundle over M x M, it has dimension n (" 2 +3 ) 
as a manifold. 

3.2. Agrachev-Sachkov formulation of rolling surfaces. A previous 
definition of a rolling map can be found in pQ, where only 2-dimensional 
manifolds imbedded into M 3 are considered. Although it only deals with the 
imbedded case, the definition of the rolling is intrinsic in the sense that it 
does not depend on the imbedding. 

The configuration space for rolling one surface on another is Q, which is 
now 5-dimensional, since M and M are 2-dimensional. A rolling is then an 
absolutely continuous curve q : [0,r] — > Q satisfying the following: if x{t) 
and x(t) are the projections of q(t) into M and M then the following two 
conditions are satisfied: 

• no slip condition: x = q(t) x(t) for almost every t £ [0, r]; 

• no twist condition: Z(t) is a parallel tangent vector field along x(t) 
if and only if q(t)Z(t) is a parallel tangent vector field along x(t). 

Notice that there is no condition corresponding to the normal no-twist, since 
the manifolds here have codimension 1. In Section H] we will show how this 
definition fits into our Definition [2j 

The no-slip and no-twist conditions can be described by means of a dis- 
tribution D in the tangent bundle of Q. By distribution, we mean a smooth 
subbundle of the tangent bundle. Then the "no slip - no twist" condition 
will correspond to the requirement q(t) £ D q M for almost every t. The 
distribution D has the following local description. In any sufficiently small 
neighborhood U C M of y £ M we pick a pair of tangent vector field e±, e 2 , 
such that {ei(x), e2(x)} is a positively oriented orthonormal basis for every 
x £ U. Define ei, e 2 in a similar way in a sufficiently small neighborhood U. 
Since the rotation group 50(2) has dimension 1, we simply need to know 
the relative angle 9 to describe q with respect to the frames given by {ei, e<i\ 
and {ei,^}. More precisely, 9 is defined by 

qe\ = cos^ei + sin6'e 2 , 
qti = — sm9ii + cos9e2- 

Thus, if Ti : Q — > M x M is the natural projection, then any q £ 7r _1 (f7 x U) 
is uniquely determined by the coordinates ( ) eU xU. 
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Let Cx, c 2 , ci and c 2 be the so-called "structural constants", defined by the 
commutation relations 

[ei, e 2 ] = ciei + c 2 e 2 , [ei, e 2 ] = ciei + c 2 e 2 . 

Define the vector fields X\ and X 2 on ix~ l (U x [/) by 

ei + cos #ei + sin #e 2 + (— ci + c± cos # + c 2 sin 6 1 ) — , 

do 

d 

e 2 — sin 6e\ + cos #e 2 + (— c 2 — ci sin 9 + c 2 cos #) — . 

Then D\ it _ 1{Ux q ) is spanned by Xi,X 2 . 

The connectivity by a curve tangent to the distribution D is the princi- 
pal problem. More precisely, given two different states g , 1\ £ <3> we as k 
whether there exists a rolling motion g : [0, r] — > Q, such that q(0) = go and 
g(r) = gi? The advantage of the formulation of no slipping and no twisting 
conditions in terms of a distribution, is that the question of connectivity 
may be reformulated through admissible sets or orbits in control theory. 

Given a distribution D on an arbitrary manifold Q, a curve q : [0, r] — > Q 
is said to be horizontal (or admissible) with respect to D if q is an abso- 
lutely continuous curve satisfying q(t) G D for almost every t. The orbit 
of D at a point go is the set of all points gi G Q so that there exists a 
curve g : [0, r] — > Q, with g(0) = go and g(r) = q\, which is horizontal with 
respect to D. We denote this set by O qo (D). It is clear that if gi G O qo (D), 
then (9 go (-D) = O qi (D). The Orbit Theorem [61 [13] asserts that C go is an 
immersed submanifold of Q and describes the tangent space of the orbit in 
terms of the diffeomorphisms of Q. A precise statement using the chrono- 
logical exponential and a broad discussion about the Orbit Theorem can be 
found in Chapter 5 of [I]. 

Also, define the flag associated to the distribution D inductively by 

D 1 = D and D t+1 = D + [D, D% 

We say that D has step k > 2 at q if k is the maximal integer, so that 

D k-i c D k = D k+i 

q q q 

If = D q for any integer k, we say that D has step 1 at q. The Orbit 
Theorem then tells us that D k q C T q O qo (D), where k is the step at q G 
O qo (D). In particular, if Q is connected and there is an integer k such 
that D k = TQ, then O qo (D) = Q. The previous result is known as the 
Chow-Rashevskh theorem [51 |9] and the distribution D is called bracket 
generating. 

We will use the expression that D has step if D has step k for any q G Q. 
Remark that if D is of step k, and there is a local basis of vector fields of 



X 1 = 

(3) 

X 2 = 
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D k in a neighborhood around any point in Q, then 

D k q =T q O qo {D). 

We now go back to the intrinsic definition presented in [I], where Q is 
the described 5-dimensional configuration space and D is spanned locally 
by ©• This definition can be restated as following: a curve q(t) : [0, r] — > Q 
is the rolling map of M on M if it is tangent to D. The main result of pQ, 
is the following description of orbits of D. Let x(x) and Se(x) denote the 
Gaussian curvature of M at x and of M at x, respectively. 

Theorem 1. For any q G Q, the orbit at qo satisfies dim O qo (D) = 2 if 
and only if x(pr M g) — x(pr^g) = for every q G O qo (D). Otherwise, 
dim O qo (D) = 5. 

Remark 4. In contrast to the definition in [Tl], the definition in p] deals 
with absolutely continuous curves. The advantage of this, is the ability to 
apply the Orbit theorem and the Chow-Rashevskii theorem. This was one of 
the reasons for us to define a rolling map in terms of absolutely continuous 
curves. Remark that all these theorems also hold if we consider Lipschitz 
curves instead of absolutely continuous. Hence, we always may interchange 
"absolutely continuous" with "Lipschitz" for all considerations in the present 
paper. 

3.3. Bryant-Hsu formulation of rolling surfaces. In [2] the authors 
give an intrinsic formulation to the problem of rolling two abstract surfaces 
M and M with respect to each other. The main tool in this formulation is 
Cartan's general method of moving frames, that is, determining canonical 
forms on an appropriate SO (2)— bundle. 

Let M and M be two connected oriented Riemannian manifolds of di- 
mension 2. Consider the respective frame bundles F, F. Then, as discussed 
in Subsection 13.11 the configuration space Q for this kinematic system can 
be identified with (F x F)/ SO(2). The conditions of no twisting and no 
slipping can be understood by means of the canonical one- forms a±, a%, a 2 i 
on F and «i,q?2, S21 on F. Recall, that these forms satisfy the structure 
equations 

dai = «2i A «2, dai = a?2i A a?2, 

da,i = -021 A«i, da.2 = — S21AS1, 

<itt2i = x «! A 02, dan = X Si A «2, 

where x and x are the Gauss curvatures of M and M respectively, see [I2| 
Chapter 7]. 

The rank two distribution D over Q corresponding to the "no slip - no 
twist" conditions is the push-forward of the vector fields, solving the Pfaffian 
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equations 

(4) a\ — ci?i = a 2 — a 2 = «2i — S? 2 i = 0, 

under the natural projection it : F x F — > Q. At the points where x—Scj^O 
the distribution D is of Cartan type, that is, the distributions 

D 2 = D + [D, D] and D 3 = D 2 + [D, D 2 ] 

have rank 3 and 5 respectively, see [2]. This implies that, under the condition 
x — x 7^ 0, the distribution D is bracket generating of step 3. To see under 
which conditions D is of Cartan type, define the following one-forms over 
the product F x F 

1 f 1 1 f 

0i = -{ai-a l ), 9 2 = -(a 2 -a 2 ), 9 3 = -(a 21 - a 21 ), 

I, ~ , 1, 
= -{ai + ai), uj 2 = -{a 2 + a 2 ), 

and observe that the following identities hold: 
6B\ = 9 3 A u 2 + -(a2i + a 21 ) A 6*2, 

d# 2 = -03 Awi - -(a 2 i + S21) A 61, 

11/ 
d9 3 = -(x — x)cji AW2 + - y(x + x)(wi A 9 2 — u 2 A 61) + (x — x)#i A 02 

Denote by £> = ker #1 D ker 6* 2 D ker 9 3 the space of solutions of the system 
OH) and let X = (X l7 X 2 ),F = (Y^Y^Z = (Z 1 ,Z 2 ) be a local basis of V 
chosen such that 

a 1 (X 1 ) = l, a 2 (X 1 ) = 0, ai(X 2 ) = l, S 2 (X 2 ) = 0, 
£*i(Yi) = 0, a 2 (li) = l > ai(y 2 ) = 0, S 2 (Y 2 ) = 1, 
ai(Zi) = 0, a 2 (Z 1 ) — 0, q?i(Z 2 )=0, S 2 (^ 2 ) = 0. 

Observe that for a sufficiently small open neighborhood U x U C M x M 
of the differential of the projection n is 

d ({ P ,c),(p,c)) n '■ T p U x 50 ( 2 ) x T pU x 50 (2) -> TpC/" x Tpf? x 50 (2) 

for any C,C E SO(2) and where T c SO(2), T e SO(2) and T c g_ 1 SO(2) are 
identified with so(2) in the usual manner. By the construction of the canon- 
ical forms on the frame bundles, it is clear that X, Y ker dir, whereas it 
is possible to choose locally Z such that Z e ker dn. Thus since ker dir has 
dimension one, we have locally 

kero?7r = span{Z}. 
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This implies that a local description of D is given by 

D = span{dn(X),dTT{Y)}. 

Recall Cartan's formula for a differential one form rj and any two local 
vector fields v, w, given by 

dr](v , w) = v(r](w)) — w(r](v)) — T)([v, w}). 

In our case, the previous equation implies the following equalities 



d9t(X,Y) = 


-di([X,Y}) 


= 0, 


d6i(X,Z) = 


-Oi([X,Z]) 


= 0, 


Mi(Y,Z) = 


-Oi([Y,Z])- 


= 0, 


d9 2 (X,Y) = 




= 0, 


d9 2 (X,Z) = 


-9 2 ([X,Z}) 


= 0, 


d9 2 (Y,Z) = 


-02([Y,Z])- 


= 0, 


d9 3 (X,Y) = 


-Os([X,Y}) 


= -fx — 5c 
2 V 


d9 3 (X,Z) = 


-e 3 ([x,z}) 


= 0, 


d9 3 (Y,Z) = 




= 0, 



It follows from these equations, that [X, Z], [Y, Z] belong to T> and [X, Y] ^ 
T> if and only if the difference of curvatures k—5c does not vanish identically. 
In fact, counting dimensions, we see that span{X, Y, Z, [X, Y]} = ker^ R 
ker 9 2 . It is clear from the choice of Z that [X, Y] ^ ker dn since if [X, Y] = 
kZ for some fcel, then d9 3 (X, Y) = —k9 3 (Z) = which contradicts our 
assumption. This implies that span{<i7r(X), dn(Y), dir([X, Y])} = Di is a 
distribution of rank 3. Analogously we obtain 

d9 1 ([X,Y],X) = -9x([[X,Y\,X]) = 0, 
d9 1 ([X,Y],Y) = -9 1 ([[X,Y],Y]) = 9 3 ([X,Y]), 
d9 2 ([X,Y],X) = -9 2 ([[X,Y],X}) = -9 3 ([X,Y}), 
d9 2 ([X,Y},Y) = 0. 
By similar considerations, we can see that 

span{X, Y, Z, [X, Y], [[X, Y],X], [[X, Y],Y}} = T(F x F), 
which implies that 

span{rf7r(X), dn(Y), dir{[X, Y]), dir{[[X, Y],X]), dir{[[X, Y],Y})} = D 2 , 

is a distribution of rank 5. 

These calculations imply that D is of Cartan type whenever x —5c does 
not vanish identically. Since the configuration space Q is 5-dimensional, 
the distribution D is bracket generating and thus, by the Chow-Rashevskh 
theorem we can completely solve the connectivity problem. In the case 
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when x = x, the distribution D is integrable and therefore Q is foliated by 
sub manifolds of dimension 2. 

It is mentioned in [2], that their construction does not depend on imbed- 
ding into Euclidean space, however no attempts are made to compare this 
definition to the one for imbedded manifolds. 

We present a simple example, illustrating the above mentioned approach. 

Example 3. Let us consider the problem of the two dimensional sphere S 2 
rolling over the Euclidean plane M 2 . We can embed these surfaces in the 
three dimensional Euclidean space M 3 via the parameterizations 

S 2 = {(cos9 cos cp, sm9 cos (p, sin cp) : — n < 9 < ir, —— < <p < —}, 

R 2 = {(x,y,0):x,yER}. 
It follows from straightforward computations that, in this case, we have 

«i = cos ipd9, a2 = dip, 0:21 = sin (pd6; 

«! = dx, S2 = dy, C(2i = 0. 

Thus, equations (jlj) take the form 

cos ipd9 — dx = dip — dy = sin >pd9 = 0. 

It is easy to see that 

da,2\ = cos ipd9 A dip = a± A «2, da 2 i = 0, 

from which it follows that x = 1 and x = 0. Since the difference of the 
Gaussian curvatures does not vanish identically, we obtain the well-known 
result that it is always possible to achieve any configuration from a given 
one by rolling the sphere over the plane without slipping or twisting. 

4. Intrinsic rolling 

4.1. Reformulation of the rolling motion in terms of bundles. Both 
formulations of rolling maps given in [T] and [2] only use the configuration 

space as a manifold of isometries of tangent spaces of M and M , without 
taking into account the imbedding into an ambient space. However, neither 
of these descriptions attempts to give any justifications for why the ambient 
space may be ignored, nor do they attempt to compare the intrinsic definition 
and the extrinsic definition given for imbedded manifolds in [TTj . We would 
like to find a reformulation of Definition [2] in such a way that the conditions 
(i')-(vi') are stated both in terms of intrinsic conditions given on Q and some 
additional conditions given on another bundle, that carries the information 
on imbedding. 

The conditions imposed over a rolling (x, g) by Definitions 1 and 2 are non- 
trivial in normal directions for the imbedding of the manifolds with codimen- 
sion v greater than 1. So, it is natural to suppose that the total configuration 
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space of the rolling dynamics will have a normal component which will takes 
care of the action of g on the normal bundle. Therefore, we make the follow- 
ing analogue construction, as we did for Q, in order to construct a principal 
bundle over MxMof isometries of the normal tangent space. We start 
from a pair of imbeddings i : M — > W n+U and T: M — > W a+I/ , given as initial 
data. Let <3> be the principal SO (u)— bundle over M, such that the fiber 
over a point x G M consists of all positively oriented orthonormal frames 
{e\(x)} v x=1 spanning T X M^ . Let $ be the principal SO (v)— bundle similarly 
defined on M. Likewise we did in Section I3~T1 identifying (Fx F)/ SO(n) 
with 

(5) Q = {qe Isom+ [T X M, T X M) x G M, x G M j , 
we identify ($ x $)/SO(z/) with 

(6) P ht :={pe Isom+ (T X M ± , T £ M X ) | x G M, x G m| . 

The space P L t is a principal SO (v)— bundle over M x M. The dimension 
of P t p is 2n + v ( v ~ x ) . The left and right actions of SO(z/) on the fibers are 

defined by the corresponding actions of SO(V) on $ and $, respectively, in 
a similar way to the left and right action on Q described in Section 13.11 We 
notice and reflect it in notations that Q is invariant of imbeddings, while 
P L p is not. 

Proposition 1. If a curve (x,g) : [0,r] — > M x Isom + (lR n+ ^) satisfies (i')- 
( vi ') of Definition [H then the mapping 

t ^ (dg(t)\ Tx(t)M ,dg(t)\ Tx(t)M ±) =: (q(t),p(t)) , 

defines a curve in Q © P h ? with the following properties: 

(I) no slip condition: x(t) = q(t)x(t) for almost every t. 
(II) no twist condition (tangential part): q(t)^Z(t) = ^q(t)Z(t) for any 

tangent vector field Z(t) along x(t) and almost every t. 
(Ill) no twist condition (normal part): p(t) (t) = ^rp(t)^(i) for any 
normal vector field ^(t) along x(t) and almost every t. 

Conversely, if (q,p) : [0, r] — > Q © P L $ is an absolutely continuous curve 
satisfying (I)-(III), then there exists a unique rolling (x,g) : [0, r] — > M x 
Isom + (M ?1+v ), such that dg(t)\ Tx{t)M = q(t) and dg(t)\ T<t)M ± = pit). 

Proof. Assume that (x,g) : [0,r] — > M x Isom + (M n+iy ) is a rolling map 
satisfying (i')-(vi'). The statements (i') and (ii') assure that 

dg(t) \t xW m e Isom(T a .(i)M, T x{t) M) and 

d<?(t)| w ^ e Isom(T ;c(t) M ± ,% ) M- 1 ). 
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Since dg(t) must be orientation preserving in M. n+U we conclude that both of 
the mappings (j7|) are either orientation reversing or orientation preserving. 
The additional requirement (vi') implies that (q,p) is orientation preserving. 
The conditions (I)- (III) correspond to the conditions (iii')-(v'). 

Conversely, if we have a curve (q(t),p(t)) in Q © P b ^ with projection 

(x(t),x(t)) into the product manifold M x M, then we may construct the 
isomorphism g G Isom + (lR n+ ^) in the following way. We write g{t) : x t— > 
A(t)x + f(i),A(t) G SO(n + u), where A(t) = dg{t) is determined by the 
conditions 

dg{t)\ Tx(t) M = q(t)\ Tx(t) M, dg(t)\ Tx(t)M ± = p(t)\ Tx(t)M ±. 

Then 

Image dg{t)\ Tx{t)M = T m M, Image dg(t)\ Tx(t)M ± = T^ t )M ± . 

The vector f(t) is determined by f(t) = x(t) — A(t)x(t). □ 

The one-to-one correspondence between rolling maps and absolutely con- 
tinuous curves in Q ©P t) f, satisfying (I)-(III), naturally leads to a definition 
of a rolling map in terms of these bundles. 

Definition 3. A rolling of M on M without slipping or twisting is an abso- 
lutely continuous curve (q,p) : [0,r] — > Q © P L> ? such that (q(t),p(t)) satisfies 

(I) no slip condition: x(t) = q(t)x(t) for almost every t, 
(II) no twist condition (tangential part): q(t)j^Z(t) = ^q(t)Z(t) for any 

tangent vector field Z(t) along x(t) and almost every t, 
(III) no twist condition (normal part) : p(t)^-^f(t) = ^-p(t)^f(t) for any 
normal vector field ^(t) along x(t) and almost every t. 

A purely intrinsic definition of a rolling is deduced from Definition [3J 
by restricting it to the bundle Q. This concept naturally generalizes the 
definition given in [1] for 2-dimensional Riemannian manifolds imbedded 
into M. 3 and we use the term intrinsic rolling for this object. 

Definition 4. An intrinsic rolling of two n-dimensional oriented Riemann- 
ian manifolds M on M without slipping or twisting is an absolutely con- 
tinuous curve q : [0, r] — > Q, satisfying the following conditions: if x(t) = 
pr M q(t) and x(t) =pr^g(t), then 
(F) no slip condition: x(t) = q(t)x(t) for almost all t, 

(IF) no twist condition: Z(t) is a parallel tangent vector field along x(t), if 
and only if q(t)Z(t) is parallel along x(t) for almost all t. 

4.2. Rolling versus intrinsic rolling along the same curves. Suppose 
that the projection of a rolling map into M x M is a fixed pair of curves. 
Questions that naturally arise are: 
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• If (qi(t) , pi(t)) and {q2(t),p2(t)) are two rollings of M on M, along 
x(t) and x(t), how do they relate to one another? How many of the 
properties of the rolling are fixed by choosing paths? 

• Suppose that an intrinsic rolling q(t) and two imbeddings, i : M — > 
~R n+u and T: M — > R n+I/ , are given. When can the intrinsic rolling 
q(t) be extended to a rolling (q(t),p(t))7 Is this extension unique? 

Before we start working with this, let us consider the following simple ex- 
ample, where the different imbeddings are easy to picture. 

Example 4. Let us consider M — R, with the usual Euclidean structure, and 
M = S 1 , with the subspace metric, when considered as the unit circle in R 2 , 
with positive orientation counterclockwise. Let x : [0, r] — > S 1 be written as 
x(t) = e tlf ^\ if : [0, r] — > R being an absolutely continuous function. Since 
SO(l) is just the trivial group, Q = M x M. It is clear from the no-slipping 
condition that 

x(t) =x{0) + <p(t)-<p(0). 

Without loss of of generality, we may assume x(0) = cp(0) = 0. We consider 
the possible rollings under different imbeddings. In the following cases, e\ 
and ei will always be positively oriented unit basis vectors for TM and TM 
respectively (when they are seen as sub-bundles of TR 1+l/ restricted to either 
M or M), while {e\} u x=1 and {e K } u K=1 are positively oriented bases of TM 1 - 
and TM- 1 . The coordinates of R 1+1/ will be denoted by (x\, . . . , x n ). 
Case 1: Let us consider the simplest example, with 

i\ : M ->■ R 2 , i x : e itp h-> (sin tp,l — cos tp), 
t i: M->R 2 , ti (z,0). 

Then 




^ = W^ 1 ^' ^ = dx~ 2 ^ 1 ^' 
Here, also SO(z/) is trivial, so there is so there is only one way to roll. 
Case 2: We do the same imbeddings as above, only increasing the codimen- 
sion by one. 

t 2 :M^R 3 , i2 '. e ilp i— >■ (sin <p,l — cos </?, 0) , 
r 2 :M-^R 3 , r 2 : x ^ (x,0,0). 
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Then 

ei (0 = cos^— (i 2 (e^)) + sin^— (t 2 (e^), 
ei (e^) = -sin^ ^-MO) +cos^ A( t2 ( e ^)), 

9 9 9 

ei(x) = — (T 2 (x)), ei(z) = — (T 2 (x)), e 2 (x) = ^— (^(x)). 

OTi OX2 OX3 

Now we know that the matrix representation B of p(t) with respect 
to the bases {e\} u x=l and {e K }^ =1 , can be represented as 

n _( (ei,p(t)ei) (ei,p(t)e 2 ) \ _ / cos0(t) sin0(t) \ , . 

V (e 2 ,p(t) ei ) (e 2 ,p(t)e 2 ) J" v -sin 0(f) cos0(t) J G bU ^' 

We calculate the restrictions of 0(i) given by (III). 

= -e{t)(smO{t)e x + cos0(f)ei) + cos 0(f) v£ (t) ei - sin 0(f) Vi (t) e 2 

= -0(f)(sin0(f)ei + cos0(t)e 2 ), 
for almost every i, so 0(i) is a constant. 
Case 5: We continue with v — 2, but change the imbedding of M to a spiral. 

6 2 :M-»1R 3 , ta : e iip h-> (sinyj, 1 - cosy?,0), 

r 3 : M — > M 3 , t 3 : x \-> —=(cosx, sinx, x). 

a/2 

Then 

ei (e"0 = cos^J-( i2 (e^))+sin^J-( i2 (e^), 
axi 0x2 

ei(e^) = -sin^ J-( t2 (e^)) +cos^ J-( i2 (e^)), 
e 2 (0 = AfcCe*)), 
er(x) = i= (-sinxj-(r 2 (x)) +cosx J-(t 3 (x)) + J^ 3 (x)) 

e x (x) = -L (-sinxJ-(T 2 (x)) + cosx A(T 3 (x)) - ^(x)) 

d d 
= -cosx—- (T 2 (xj) -sinx-—(? 3 (x)). 

OX\ OX2 
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We have the same matrix representation of p(t), 

n _( (ei,p(t)ex) (e 1 ,p(t)e 2 } \ _ f cos6(t) sin0(t) \ , - 

(e 2 ,p(t) ei ) (e 2 ,p(t)e 2 ) )'{ -sm9(t) cos0(t) J G bU ^' 

We calculate the restrictions of 9(t) given by (III). 

D L 

P(t)Vi {t) ei = = -jjrPitfa 

= -0(t)(sin0(t)ei + cosS(t)ei) + cos0(f)VjLei - sin0(t)V± (t) e 2 
x(t) 

so 9(t) = 9q + ^=sc(t). So now, the circle M will rotate along the 

spiral M, but its path is determined by the initial angle. Notice also 
that if we define a new orthonormal frame of TM L by 

Ti = cos I -= J e x - sin ( -= j e 2 , 



6(t) (sin 0(t) e a + cos0(t)e 2 ; 



\/2/ V\/2/ 

T 2 = sin(^)e 1 + cos(^)e 2 , 

then becomes a constant matrix with respect to the bases e±, e 2 
and Ti, T 2 . 

We see that for cases above, the intrinsic rolling t i— > (e lip ^\ (f(t)) either 
uniquely induces a rolling, or the rolling is determined by an initial config- 
uration of the normal tangent spaces given by 0(0) = 9q. Note also that we 
are able to find a choice of bases so that p(t) is constant with respect to this 
basis. Notice that these bases consist of normal parallel vector fields. 

We continue to work with oriented manifolds M and M imbedded in M. n+1/ 
and containing curves x(t) and x(t), respectively. In the remaining of this 
section we will use the following notations: {ej(t)}™ =1 will be a collection 
of parallel tangent vector fields along x(t) that forms an orthonormal basis 
for T x < t \M at each point of M, {e\(t)}^ =1 will be a collection of normal 
parallel vector fields along x(t) forming an orthonormal basis for T x < t \M . 
We know that we can construct such vector fields by parallel transport and 
normal parallel transport along x{t). Parallel frames {ej}" =1 and {e K } v K=1 
will be defined similarly along x(t). Recall that Latin indices . . . always 
go from 1 to n, while Greek ones k, A, . . . vary from 1 to v. 

The following lemma reflects that a rolling map preserves parallel vector 
fields. Namely, the image of a parallel frame over M has constant coordinates 
in a parallel frame over M. 
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Lemma 2. A curve (q(t),p(t)) in Q © P L $ in the fibers over x(t)), 
satisfies (II) and (III) if and only if the matrices A(t) = (aij(t)) and B(t) = 
(b KX (t)), defined by 

a i3 (t) = e*(t)g(t) ej (t), b KX (t) = e* K {t)p(t)e x (t), 

are constant. 

Proof. Let (q(t),p(t)) be an absolutely continuous curve. Then we have 
(ei,ej) = (ei,ej) = and 



by the product rule. The vectors q(t) 1 ej, q{t)ej are tangent, so (q(t) 1 e i , ef) = 
(e-i, q(t)ej) = and 

aij(t) = (ei,q(t)ej) + (ei,q(t)ej) + (ei,q(t)ej) 

ei >dt q{t)e i- q{t) dt ej )=°- 

So (II) holds if and only if a^ (t) = 0. Similar result holds for the basis of 
the normal tangent bundle. □ 

The following two theorems give the answer to the questions raised at the 
beginning of this section. 

Theorem 2. Let qo '■ [0, r] — > Q be a given intrinsic rolling map without 
slipping or twisting with the projection pr Mx ^q (t) = (x(t),x(t)) . Denote 
by k the dimension of the space of parallel tangent vector fields along x(t), 
that are orthogonal to x(t) whenever they are defined. Then the following 
statements hold. 

(a) The map qo is the unique intrinsic rolling of M on M along x(t) and 
x(t) if and only if k < 1 . 

(b) If k > 2, then there exists an infective Lie group homomorphism 
C : SO(fc) -> SO(n) such that for each A' e SO(fc) the map q (t)((A') 
is an intrinsic rolling over (x(t),x(t)) , and any intrinsic rolling over 
(x(t),x(t)) is of this form. 

Proof. Pick up frames of parallel vector fields {ej}" =1 and {ej}™ =1 along 
x(t) and x(t), respectively, such that qo(t)ei = e^. This is possible due 
to Lemma [2j We also choose the frames in a way that the k first vector 
fields are orthogonal to x. 
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Writing x = YJi^iif)^) and * = Ya=i ^(X) e i(X)> we g et ^(i) = ±i(t) 
and x\{t) = ■ ■ ■ = Xk(t) = 0. So, if q is any other rolling, then A = (a^) = 
((e,(t), q(t) ej(t))) is clearly of the form 



where l n -k is the ((n — k) x (n — A;) ) -unit matrix. 

The converse also holds; that is, for any such matrix A, there is a rolling 
corresponding to it. □ 

Theorem 3. Let q : [0, r] — > Q be an intrinsic rolling and let i : M — > 

M n+!/ and T : M — > M. n+l/ be given imbeddings. Then there exists a rolling 
(qo,p) '■ [0,t] — > Q © P h ? that is unique up to a right action of SO(u). 

Proof. We pick up normal parallel frames {e\(t)}^ =1 and {e K {t)} v K=1 along 
x{t) and x(t), respectively. For any element B G SO(i^) we define pit) by 

B = (e K ,p(t) e x ). 

The map (q (t),p(t)) is a rolling by Lemma [21 and all rollings are of this 
form. □ 

Corollary 1. Assume that x(t) is a geodesic in M. Then there exists an 
intrinsic rolling of M on M along (x(t),x(t)) if and only ifx{t) is a geodesic 
with the same speed as x(t). Moreover, if n > 2 then there is an injective 
Lie group homomorphism (: SO(n — 1) — > SO(n), such that all intrinsic 
rollings over (x(t),x(t)) differ by an element in Image (. 

Proof. Taking into account the equality ^x(t) = ^q(t)x(t) = q(t)^x(t), we 
conclude that if x(t) is a geodesic then x(t) is also geodesic. In order to 
satisfy (I) we need to require that the speed of x(t) is the same as the speed 
of x(t). Conversely, the equality of speeds implies condition (I). 

We start the construction of rolling map by choosing e\ (t) = t^t^t^ that 
is parallel along x(t). The remaining n — 1 parallel vector fields we pick up 
in a way that they form an orthonormal basis together with e\(t) along the 
curve x(t). We repeat the same construction for a parallel frame {ej(t)}" =1 
along x(t). Define the intrinsic rolling q(t) by 

e\(t)q(t) ej (t) = e*(t)q(t) e x {t) = 5 ld , 

(8) A'^ie^q^e^it))^, 

where A' £ SO(n— 1) will be a constant matrix. Conversely, we can construct 
a rolling by formulas flH]) starting from A' £ SO(n — 1). □ 



24 



M. GODOY, E. GRONG, I. MARKINA, F. SILVA 



5. Distributions for rolling and intrinsic rolling maps 

The aim of this Section is to formulate the kinematic conditions of no- 
slipping and no-twisting in terms of a distribution. In this setting, a rolling 
will be an absolutely continuous curve almost everywhere tangent to this 
distribution. 

5.1. Local trivializations of Q. Let 7r : Q © P b $ — >■ M x M denote the 
canonical projection. Consider a rolling j(t) = (q(i),p(t)), then n o = 
(x(t),x(t)). Given an arbitrary t in the domain of 7(t), let U and U denote 
neighborhoods of x(t ) and x(t ) in M and M, respectively, such that the 
both bundles TM — > M and TM X — > M trivialize being restricted to U. 
In the same way we chose U, such that both TM — > M and TM 1 - — > 
M trivialize when they are restricted to U . This implies that the bundle 
Q®P l ,t — > M x M, trivializes when it is restricted to U x U . To see this, let 
{ej}™ =1 , {e\}^ = i, {ii}^ =1 and {e K } y K=1 denote positively oriented orthonormal 

bases of vector fields of TM\u, TM- 1 ^, TM\q and TM ± \jj, respectively. 
Then there is a trivialization 

(9) Q®P^\ Ux o A i7x(7xSO(n)xSO(i/) 

(q,p) h-> (x,x,A,B), 

given by projections 

x = pT u (q,p), x = pr (q,p), 

A = ( a n)lj=i = ((^i, e l ))'l j=1 , B = (6« a )«,a=i = ((P £ A, e K ))«, A =i • 

The domain of 7 can be chosen connected, containing t , and such that its 
image lies in 7r _1 (?7 x U). Let us identify 7(2) with its image under the 
trivialization given by (x(t),x(t), A(t), B(t)). 

Each of the requirements (I)- (III) can be written as restrictions to j(t). 
We will show, that all admissible values of j(t) form a distribution; that is 
a smooth sub-bundle, of T(Q © P h i}- We will use the local trivializations to 
describe this distribution. 

5.2. The tangent space of SO(n). Let U and U be as in Section I5TT1 Then 
we get in trivialization 

Tn-\U xU)=TUxTUx TSO(n) x TSO(u). 

The decomposition requires that we present a detailed description of the 
tangent space of SO(n) in terms of left and right invariant vector fields. 

We start by considering the imbedding of SO(n) in GL(n), the group of 
invertible real n x n matrices. Denote the matrix entries of a matrix A by 
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(djj) and the transpose matrix by A 1 . Then, differentiating the condition 
A 1 A = 1, we obtain 

n 

T SO(n) = P| ker w^-, = (a r j da ri + a ri da n ) . 

i<j r=l 

It is clear that the tangent space at the identity 1 of SO(n) is spanned by 

d d 

W i:j (l):=- — , l<i<j<n. 

Odij Odji 

We denote so(n) = span{Wij(l)} following the classical notation. We use 
the left translation of these vector to define 



(10) 



n ( d d \ 



as global left invariant basis of T SO(n). Note that the left and right action 
in TSO(ra) is described by 

— - V — JL a-Y" — 

dan " da r j dan , ^ da is 

j r=l J J s=l 



We have the following formula to switch from left to right translation 

A ■ — — y a r i~ — / a riO],l-^ — / O-riO'siO'slT: 

octii oa r j oa r i oa r i 

%3 r=l r J l, r =l " ,r,s=l 



and the other way around, 



T, ' A — y Ojj s — — y CLjsOi l— — y Oj'sOir-O/rT; — 

dan , da is ^ oa [s ^— ' <9a* s 

»=1 I,s=l «,r,s=l 



Therefore, the right invariant basis of TSO(n) can be written as 
W y (l) • A = Ad(A" 1 )l^ i (A) = ^(a ir a is - a ir ai S )W rs (A). 



If we let Wij be defined (fTUjl also when i is not less then j, (so Wij = —Wji) 
then the bracket relations are given by 

[W^, W kl ] = S jtk W a + 8 itl W jk - S itk Wji - 5 jtl W ik . 
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5.3. Distributions. Now we are ready to rewrite the kinematic conditions 
(I)- (III) as a distribution. Let 7(2) be a rolling satisfying the conditions 

(I) - (III). Consider it image under the trivializations. Then 

n d v d 

(II) T (t) = ±{t) + s(t) + a l3 — + h * x db^- 

i,j=l l 3 ft,A=l 

If we denote x(t) by Z(t), then (I) holds if and only if x(t) = q(t)Z(t). 

We want, basing on conditions (II) and (HI), write the last two terms 
in f TTTT) in right invariant basis of corresponding tangent spaces of SO(ra) and 
SO(z/). We start from (II) and remark that 

n n 

<?(X) e j = ^2 a ij(t)eu and q~ l (t)ei = S ^a ij (t)e. j 
i=i j=i 

for orthonormal bases {ej}^ =1 and {e}}" =1 . Condition (II) holds if and only 
if q%ej(x(t)) = %qej(x(t)) for j — 1, . . . , n, that yields 

n n 

= ^2 an (Vz(t)ej, e;) - ay - ^ ay (V g z(t)ej, e f ) 

for every z,j = 1, . . . ,n. Hence, the third term in ( II ip can be written as 
follows 

n q n / n n \ d 

i,j=l l i i,j=l \i=i 1=1 ) % i 



n d n d 

J,i=l i,Z=l 



l,J = l J l,],r,S=Y 



dan . J " '^-"-v" ^ 

11 1 ft j fi \ \ q 

i,j=l \ s=l \ r=l I J % 3 

a 



(12) = « V ^W e i' e i) ~ ( v ff*(*)? e i> <?e*)) A 



ij=i J 
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The coefficients in the basis A - in the sum (fl~2j) are skew symmetric, from 



the property of the Levi-Civita connection. Now we can write 

_d 

da 



n d 

(13) Yl h ^-^r. = Yl (( v ^w e i' 6i ) ~ ( v iz{t) ae ^ ae i)) w ai A )- 



i,j=l l 3 i<j 



Written in a right invariant basis, we obtain 
(14) 

n d 

d v]fa7. = « v ^)<r 1 e J) g~ 1 e i ) - (Vgzftjej, ei» Ad^" 1 )^-^). 

i,j=l l ^ i<j 

Similarly, (III) holds if and only if 
(15) 

E kxst x = E (( V z ( t) e A,e K ) - (v^ (t) pe A ,pe K \) W kX (B). 

K,X=1 K<\ X X ' N ' 7 

E ((V^-^A.P- 1 ^) - (vfe (t) e A ,e K )) Ad( J B- 1 )W KA ( J B). 



Definition 5. //X zs a vector field on M , then let us define V{X) and 
V (X) the vector fields on Q © P L x, such that under any local trivialization 
h as in flU]) and any (q,p) E 7r _1 (x) they satisfy 



(16) dh{V(X){q,p)) = ((Vx W e 3l e 4 ) - (V qX{xm , qa}) W tJ {A). 

i<j 

(17) dh (V^X^p)) = £ «Vi (x) e A ,e K > - (V^ x{x) pe x ,pe K )) W kX {B). 

Notice that if Y(x) = X(x) = X E T X M, then V(Y)(q,p) = V(X)(q,p) 
for every (q,p) E (Q © P h t) x - Hence, we may define V(X )(q,p) whenever 
X E T x M and (q,p) E (Q © P h ?) x . Also notice that the map X H> V(X) is 
linear. The same holds for V x . 

Remark 5. Notice that, at first glace, it may seem that all of the coefficients 
of Wij(A) and W K \(A) in (ITBl) and (1T7|) vanish from conditions (II) and (III). 
This is not true, however. Even though, for any tangential vector field X 

jX{x{t)) = V m X(x(t)), 

in general, V qx ( t )q(t)ej does not coincide with ^q(t)ej(x(t)). To see this, 
notice that 

f t a sj e s (m) = a s Mm)+a SJ V m eMt)) = d,e s (x(t))+a, V„ (i) e s (x(t)), 
while 

Vg i(4 )a si e s (x(t)) = a sj V qx{t) e s (x(t)). 
Similar relations hold for 

at 
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We may now sum up our considerations that have been made in this 
Section in the following result. 

Proposition 2. A curve (q(t),p(t)) in Q © P L $ is a rolling if and only if it 
is a horizontal curve with respect to the distribution E, defined by 



E {qjP) = {X + qX + V(X )(q,p) + V ± (X )(q,p)\ X G T X M) , (q,p) G (Q©P t ,r) 



If we use the same symbol to denote the restriction of V(X) to Q, we also 
have 

Proposition 3. A curve q(t) in Q is an intrinsic rolling if and only if it is 
a horizontal curve with respect to the distribution D, defined by 



5.4. Properties of the distribution. We present some of the properties 
for the distribution E that basically reflects the results found in Theorem [2] 
and El 

Proposition 4. (a) E is biinvariant under the action of SO(^). 

(b) Let X be a vector field on M. If for any q G Q x ,x, if (A q)X(x) = 
qX(x), then 



Proof. To prove (a) , we only need to show that B - V 1 - (X) (g, p) = V 1 - (X) (q, B p) 
and V L {X){q,p) ■ B = V ± {X){q,pB ) for any B = (M^i e. First, let 
h be a local trivialization as in ([9]) and (q,p) G (Q © P L ^) X - Then, from ( |T5l) . 
to show V ± (X)(q, B p) = B ■ V ± (X)(q,p), it is sufficient to show that if 



D q = {X + qX + V(X )(q)\X G T X M} , q G Q x . 



V(X){A q)=A V(X)(q). 
Similarly, if (qA )X(x) = qX(x), then 

V(X)(qA )=V(X)(q)A . 




then 



V 




Since 




V 



V 



( 



U V 



we need to show that 
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To obtain this, let Y be any exstension of the vector field qX, to M, n+U . Let 
e M and e# be any exstension of e M and e#, and write 

n+i> q 

^ = £/*(*) si? 
i=i 

Then 

ia=l i9=l 

= (VyBoT^Bo^) = (^o^F/i^)— ,B&> \ = (BoVyB T^Bo^) 
We show right invariance by showing that 

V 

in a similar way. 

The proof of (b) is totally analogous to the proof of (a) □ 

6. A CONTROLLABLE EXAMPLE: S n ROLLING OVER R" 

6.1. Formulation of the rolling. We want to illustrate the properties of 
the distributions, by proving that the unit sphere S n in IR n+1 rolling over R n 
is a completely controllable system, by showing that the distribution D is 
bracket generating. This result was obtained in [TJ] , but we want to present 
this example here in order to illustrate the advantages of our theory. 

Consider the unit sphere S n as the submanifold of the Euclidean space 
R n+1 , 

S n = {(x , ...,x n )e R n+1 \ xl + --- + x 2 n = l}, 
with the induced metric. 

For an arbitrary point x — (^0; • • • j ■En) £ S n , at least one of the coordi- 
nates Xq, . . . ,x n does not vanish. Without lost of generality, we may assume 
that x n 0, and consider the neighborhood 

U = {(x ,...,x n )eS n \ ±x„>0}, 

where the choice of the ± sign depends on the sign of x n . To simplify the 
notation, we define the following functions on U 

n 

S ji X ) = ^ ] x r ■ 
r=j 

These functions are always strictly positive on U, and we use them to define 
an orthonormal basis of TU . We will write simply Sj instead of Sj(x), since 
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dependence of x is clear from the context. Define the following vector fields 
on U 

These vector fields form an orthonormal basis of the tangent space over U. 
d 

We set ii = — to be the standard basis of M. n . 

Before proceeding with the necessary calculations, let us state two tech- 
nical Lemmas whose proofs can be found in section 16.21 and 16.31 

Lemma 3. Let 1 < i < j < n. Then 

(^efe e j' e i) = , = ~~ (^e fe e i; e j) ; 

for any k — 1, . . . , n. 

Remark 6. The properties of the connection V have the following conse- 
quences: 

• The compatibility of V with the metric and (e^, ej) = 5ij, imply that 

(V ek ej,ei) = - (V ek ei,ej) . 

In particular, (V efe ej, e,) = 0. 

• The symmetry of V, imply that if I < k, then 

n 

[e k , ei] = V ek ei - V e; e fc = V" (V efe e; - V e; e fc , = -1 e k . 
Lemma 4. For k, I = 1, 2, . . . , n 

k > I 

k = l 

Sk 

- , k<l. 

^/Sfc-lSfeSfcSj-l 



e k 



Xl-l 



y/Sl-li 



It is a direct consequence of the choice of the vector fields e k that Ve fc e/ = 
0, and [e k , e{\ — for all k, I — 1, . . . , n. 

Consider the vector fields X k = e k + qe k + V(e k ) which generate the 
distribution V, introduced in Proposition [31 restricted to U. In this case, 
we have the explicit form 

n k—1 

X k (x, x, A) = e k (x) + ^ a ik^i{x) - ^2 W ik (A). 

t=i i=i v Si ~ lSi 
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In order to determine the commutators [X k , Xi], let us assume that k > I. 
Then 

k— 1 n l—l n 

[X k , Xi] = [e fc , e{\ - y~] E Wikdjiij + J ~ Wjia ih ei 

i=l j=\ V * 1 4 j_i j=i V J J 



X> ( J=) -*+|> (^) ^ft-^-^ 



fc— 1 n (— 1 n 



V ' {=i 3=1 V 1 j-i j=i V -J -1 J 

-\w lt - £ w. + E E ;■ . (- w + 

(n fe-1 !-l 

e fc + £ " E ~F^=W* ~ E "F^i* 

1 x 2 

— w ik + E 



ti v^T^ ) \ s t Sl _ lSl \ Sj Sj 



V \ 3 =1 r=l V / \ 3 = 1 

Xl-1 



X k -W t 



Ik- 



y/Si-iSi 

Define the vector fields Yi k , for I < k, by 

^zfe : = [^ij^fe] H . X k = Wik- 

y/Sl-lSl 

Finally, let 

n 

Z 1 = [Y 12 ,X 2 ] = y^agej, 

i=l 

n 

Z k = [X x , Y lk ] = ^ a ik&i , k = 2,...,n. 

i=i 

We conclude that the entire tangent space is spanned by {X k }^ =1 , {Y lk \i<i <k < n 
and {Zk}k=i- Hence, D is a regular bracket generating distribution of step 
3, which implies that the system of rolling S n over W 1 is completely control- 
lable. 
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6.2. Proof of Lemma [31 The proof of this Lemma is rather technical and 
it consists mostly of rewriting formulas in an appropriate way. We begin 
with some observations. 

• s j-i = %2 j-i + s j- 

• If H is the Heaviside function 



H{x) 

then 



1 when x > 
when x < ' 



d 

-sj = 2x k H(k-j). 



dx k 

For any integer j, 

H(j) = 6 0>j + H(j-1). 
Due to the identity 

d \ I Si ( x k x^iH(k - i) 
-d k ,i-i H 



dx k I V 
we obtain 

. A e .\ _ /Ay (_ x 5 . 1 + ^-isffifr ~ 



Si-! V s « 

, "sT . s k H(k-i- 1) 



Si—! \ 



Sfep 1: Finding V _a_ej. We calculate 

d I sj x k H{k — j) , . . / Sj 



dx k Y Sj-! y/Sj-lSj 



- x k H(k - j + 1) 



3 



x k H(k- fiJ— 2 - [ — — ) - Xkfaj-iJ^ 

— / x k x)_ x H(k - j) x k 5 kJ -! 



Sj_i \ s j-l s j s j-l 

and get 

x k x\,H(k-3) x k b kA , 
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x k x]_ x H{k - j) 2x k H(k - j) x k 8 k j- X ' 



Sj-iSj Sj Sj-i 



s j (fik,j-ix^ x d 2x k H(k-j) d t x j ^ 1 H(k-j) d 



Sj — i y Sj . doc f s j doc j — \ s j doc ^ 



it/1 -\ s i-i + s i i x k8k,j-i\ 
x k H(k - j)^ J - H ^— I ej 

17^ 2x fc - 



Step 2: Calculating V efe ej. Using Step 1 and formula f lLSj) . we are able to 
compute 



'I S k -1 \ S k ^ &*! 



V s fc-i VV s j-i s j s j-i J 

-7^ ^% Xr ~k. + H{k ( x ^dxh ~ 2Xk - 1 ^)) 

S k V V Sj-iSj Sj--! J 

f iy5r j " liJ(J " A; " 1) ^ Xr ^ +Xj " 1 ^ ^- 2w{ ^ 

J ' J \ r=i r=mnrx{ i.k\ 



Step 3: Obtaining (V es .ey, e,). We calculate it case by case, 
• if = j , then 
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e fc J 2_^x, 




Sk-lSk S k V S k -1 ^—f 0x r 

r=k 



Sk d 
s k-i dxk-i 



1 " 

— E 



9 



r=k— 1 



and so 



(V efe e fc ,ei) 



■Ei—l / ^ 



if k < j, then 



<5i,fc-i + 



Sk-\H{k — i — 2) \ Xi-\H{k — i — 1) 



gfc ( f _ ( s i-l + , x j-l 



+— 



1 / Si 



d 







s j V s i~ 1 



dx r 



r=3 



r=j 



Sk I Xk-1 



Sk-1 \ Sk 

if j < k, then 



-2e J +2e J )) =0 



5fc 



s i-i + s j i „ 1 / s i 



Sj-lSj 



e i - 



s i V 



5 



dxj_i 







2x fc -i- 



Xk-l 
Sk 



Sk' 



S 3-l + S j 
Sj-lSj 



e,- H 



1 / s., 



s i V ^i -1 



r=k 



dxj^i 



I s k Sj I Xj-i d 2x k -i d Xk-iXj-i x _g 2x fc „i <9 

Sk-iSj-x \ Sj dxk-i Sj dxj-i SkSj r dx r Sj dxj^i 

Xj-i 

-e k 



y/Sj-lSj 

The conclusion is that all the Christoffel symbols T\j = (V efe ej, e;) vanish, 
except for 
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6.3. Proof of Lemma HI We need to consider three cases. Observe that 



dx k \y/Sl-iSlJ y/Sl-!Sl Sl-i y/SnSl 

1 x^iX k H(k -I) 1 x^iX k H(k - I) 



Si^i y/Si-iSi Si y/St-iSi 



1 / (si + s^ 1 )x l ^ 1 x k H(k - 1) 

0\-X,kSl 



so 



xi-i ^ _ I s k Y ^ s ( s i + si-i)xi^ 1 x k ^ 1 H(k - I - 1) 



n / 



s/c " V si 

r=k 



Sk f , . (si + s^ 1 )xi^ 1 x k ^ 1 H{k - I - 1) 



(z,_ lS ^(Z - 1 - k) - {SI + 8 ±^t^m 

If fc = Z, then 

e fc ( , = -*— -s k 

sJSk-lSkJ 4-1 V s fc s fc 

1 s\ + (s fe + - s k ) 1 



If k > L then 



\ / s k ( (si + s;_i)a:;_iXfc-i 3:^1 (g; + g;_i)a:;„is fc 



• If < Z, then 

k \y/Sl~lSl) \l sf^SiSk-l \ S k \ 1 Si 

Xk-iXi-i 

y/Sk-lSkSl-lSl' 
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7. A NON-CONTROLLABLE EXAMPLE: SE(3) ROLLING OVER M 6 

7.1. Calculation of the dimension of the orbits. Let SE(3) be the 
group of orientation preserving isometries of M 3 . We consider the case of 
SE(3), endowed with a left invariant metric defined later, rolling over its 
tangent space at the identity TiSE(3) = se(3), with metric obtained by 
restricting the left invariant metric on SE(3) to the identity. Our goal is to 
determine whether any two points in the configuration space can be joined by 
a curve tangent to the distribution presented in Definition [5j This problem 
is equivalent to the controllability of the system, that is, we want to obtain 
any configuration by rolling without twisting or slipping, from a given an 
initial configuration. 

We give SE(3) coordinates as follows. For any x G SE(3) there exist 
C = (cij) G SO(3) and r = (r!,r 2 ,r 3 ) G M 3 , such that x = (C, r) acts via 

x{y) = Cy + r, for all y G K 3 . 

The tangent space of SE(3) at x — (C,r) is spanned by the left invariant 
vector fields 

ei = Yl = 7l{ C 'dL~ C ' A) 

(19) 



e 2 = Y 2 



(20) 



_L ( r -JL 

V2 1° ' den 



c ® 

desi 



v/2 



C~i-$ 



33 



c i 3 at 



ji 



e 3 = Y 3 



(21) 



1 / £j _ d q _ d \ 

\/2 \ 9c23 dcz2 ) 



J_ v^3 
V2 2-'3= 



r 9 



■33 



J 2 



(22) e k+3 = X k = C- — = J2c 3k -^r fc = 1,2,3. 

j=l 3 

Define a left invariant metric on SE(3) by declaring the vectors e±, . . . , e§ 
to form an orthonormal basis. The mapping 
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permits to identify se(3) endowed with the induced metric, with M 6 with the 
Euclidean metric. We write = tJt on M 6 and try to see how the intrinsic 
rollings of SE(3) on M 6 behave. Note that Q = SE(3) xtfx SO (6), because 
both manifolds SE(3) and IR 6 are Lie groups, so their tangent bundles are 
trivial, and dim Q = 27. 

Let us denote by V the Levi-Civita connection on SE(3) or M 6 with re- 
spect to the corresponding Riemannian metrics defined above. The covariant 
derivatives V ei ej are nonzero only in the following cases 

V Yl F 2 = -Wy 2 Y 1 = -^=Y 3 



V Yl Y 3 = -V Ya Y 1 = - 7 =Y 2 



Vy 2 F 3 = -Vy 3 F 2 = -^=ri 



1 

1 

2\[2 



Vy 2 Xfc = —j= (S 3 ^Xi — 5i t kX 3 ) 

Vy 3 X fc = — — (63^X2 — &2,kX 3 ) , 

where 5y denotes the Kronecker symbol. On the other hand, it is well- 
known that Vg 4 ej = for any Proposition [3] and Definition [5] show that 
the distribution D over Q is spanned by 

Z 1 = Y 1 + qY 1 + ^W 23 + ^W^ 



Y 2 + gY 2 - ^W l3 + j-W*, 



(23) 



'56, 



Z 3 = Yt + qYs + ^Wu + j-W 

K x = X 1 + qX 1 , 

K 2 = X 2 + qX 2 , 

K 3 = X 3 + qX 3 . 

In order to determine the controllability of rolling SE(3) over M 6 , we 
employ the Orbit Theorem [61 [13]. In the case of D, defined by the vector 
fields fT23|) straightforward calculations yield that the flag associated to D is 
on the form 

D 2 = span {^12,^13,^23}, 
(24) D 3 = D 2 ®span{qY 1 ,qY 2 ,qY 3 }, 

D 4 = D 3 , 

and so dim/} 2 = 9, dimD h = 12 for all k > 3 and the step of D is 3. 
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Let (x ,x , A q ) be an arbitrary point in Q, and let O 



(x ,x ,A ) 



denote the 



subset of all points in Q which are connected to (x ,x , A ) by an intrinsic 
rolling. The Orbit Theorem asserts that, at each point, D 3 is contained in 
the tangent space of the orbits. However, since we know that D 3 has a local 
basis, we have the stronger result of 



T, 



(x,x,A)O(x ,x ,A ) 



n 3 



holding for all (x,x, A) e O^ ,a )- 

It follows from ( 1211) that O( X0 ^ 0j a ) has dimension 12. Since 0^ xo ,x 0> Ao) is 
not the entire Q, we conclude that the system is not controllable. 

We end this Section with a concrete example of an intrinsic rolling q[t) = 
(x(t),x{t),A(t)), where 

x(0) = id R a, x(0) = 0, A(0) = 1. 

Define the curve x : [0, r] — > SE(3) by 



(25) x(t)y 



I cos9(t) sm9(t) 
-sm6(t) cos^(t) 
\ 1 




where 9(t) and i/j(t) are absolutely continuous functions with 8(0) = ip(0) = 
0. Then x = \f29(t)Y l + ip{t)X 3 for almost every t, and the rolling has the 
form q = \[2Q{t)Z\ +ip(t)K 3 , or equivalently 



(26) 
(27) 

(28) 



x(t) 
x(t) 

A(t) 



V2e(t)Y 1 + ^p(t)x 3 

V2 6(t) qY 1 + ijj(t) qX 3 



6{t) [ -W 23 (A) + W 45 (A) 



for almost every t. It follows from equation ( 12 7p that 



x(t) = 

Equation fl28|) can be written as 
A(t) = A >W 



( V20{t) \ 







W 23 (l) + 6(t) W, 
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which can be solved by exponentiating to obtain 



A(t) 



( 1 









cos (V) 
• ( sit) ' 

o 






sm 





cos i^-Y 











cos0(t) 
-sm6(t) 









sin 6(t) 
cos 9{t) 




0\ 








1 J 



7.2. Imbedding of SE(n) into Euclidean space. Since it is less obvious 
how to extend an intrinsic rolling of SE(3) on se(3) to an extrinsic rolling 
in ambient space, we describe an isometric imbedding of SE(ra) into the 



Euclidean space 



(n+l) 5 



. Identify an element C G IR ( - n+1 ) with the matrix 



C 



Cll 



Cn+1,1 




Define the inner product on R( n+1 ) 2 by 

<Ci,C 2 >= trace ((Ci)*C 2 ). 

Note that since 



71+ 1 



the metric (•, •) coincides with the Euclidean metric. From this we get that 



n+l 



den , 

11 J 1,3=1 



is an orthonormal basis for the tangent bundle TIR*™" 1 " 1 ) 2 with 

respect to (•,•). 

We define the imbedding of SE(3) into R( n+1 ) 2 by 



/ : 



x 



SE(n) 
:(C,r) 



)(n+l) 2 



h+ C 



C r 
1 



This mapping is in fact an isometry of SE(n) onto its image. To see this, no- 
tice that the metrics coincide at the identity, and that the metric of IR( n+1 ) 2 , 
restricted to Images, is left invariant under the action of SE(n). Hence, the 
metrics on SE(n) and Image i coincide, and i defines an isometric imbedding. 

7.3. Extrinsic rolling. We will use the imbedding from Subsection 17.21 to 
construct an extrinsic rolling of SE(3) over sc(3) in M 16 . We use cV,- to 
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denote For the sake of clarity, we denote by M the image of SE(3) by 



/ Tnpri 1"np A/ppfoy fiplnc; 

u. 1 11C11 ij±±C VClilWl 11C1U.D 


spanning TM are 




ei = 


Yx = 




- c i2 da) , 


e 2 = 


Y 2 = 


3 

V 1 = 1 


- Q 3 9ii) , 


e 3 = 


Y 3 = 


1 3 

* i=l 






Xk = 


3 


fc = 1,2,3, 






3=1 





where we suppressed dt in the notation. We introduce an othonormal basis 
of TM X 

1 3 

T l = -7= Yl ^' 1 ^' 2 + C V d ^ ' 
V 7=1 



T, 



1 

1 3 



= A = 1,2, 3 

2 A = 9 4m , /i = 1,2,3,4 
We denote by M the image of M 6 into R 16 by the imbedding 



71 £l 







i 



(Xi,X2,X 3 ,X4,X 5 ,X 6 ) H> 

V o o 

We have the following orthonormal basis of TM 



72^3 X 5 
1 X 3 Xq 



/ 









ei = 




- 9 2 i 


e 2 = 


>• 


-d 31 


e 3 = 


^3 


- d 32 
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— <9/c4 

while the vector fields spanning TM 

1 

i 

71 



(2 



k = 1,2,3 

(<9i 2 + <9 2 i), 

(9i 3 + 5 31 ), 



£3 



-^=(d 23 + <9 32 ), 



e 3+K — 9 KK , K — 1, 2, 3, 
= ^4k k — 1) 2, 3, 4. 
In order to extend an intrinsic rolling q(t) with n(q(t)) = (x(t),x(t)), we 
will find an orthonormal frame of normal parallel vector fields along x(t) 
and x(t). Along x(t), we may use the restriction of {e^}^. For the curve 
x(t) the answer is more complicated. 

We first study the value of V x for different choices of vector fields. 

(1) V^Ha = 0, for any tangential vector field X. 

(2) V^ fc T = 0, for any normal vector field T. 

(3) Otherwise 





Ti 


T 2 


T 3 


*1 


*2 


^3 






i y 

2V 7 ! 3 











vi 


i >v> 

2^ 3 


-* 3 ) 


1 y 







|T 2 




I "Y 

~i72 l2 


1 


| (*2 - * 3 ) 







|T 3 



We use the above relation to construct an extrinsic rolling. We will illustrate 
this by considering the curve (|25|) . 

Since x(t) = y/2 9{t)Y 1 {x{t)) + ip(t)X 3 (x(t)), the vector field 

3 

*(*) = E M*)T A (z(t)) + t; 3+A (t)* A (x(t))) , 



A=l 



is normal parallel along if 
vx - -j= (v 4 - v 5 ) J Ti + ( i> 2 + j T 2 + I v 3 - -v 2 ) T 3 



V4 + A^Vl ) #1 



T> 5 - -^=Wi ) ^ 2 + V 6 V 3 = 0. 



Hence we define a parallel orthonormal frame along x(t) by 

1 1 

ei(t) = cos0Ti(x(t)) --=fun9Vi(x(t)) + -=sm6y 2 (x(t)) 

v 2 v 2 
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e 2 (t) = cos f T 2 (z(t)) + sin f^j T 3 (x(t)) , 
e 3 (t) = - sin [0 T 2 (x(t)) + cos f^j T 3 (x(t)) 



e 4 (t) = -L sinflT^t)) + ^Ltlv^t)) 



cos 9 



*2(x(t)), 



1 . ™ / , ,\ 1 — cos Q -r , , \ \ 1 + cos 6 1 T . ... 
e 5 (t) = --j= sinflT^t)) + *i(x(t)) + * 2 (x(t)) , 

e 6 (*) = tt 3 (s(*)) , 
e 6+A (t) = H A (x(t)) , A = 1,2,3,4. 

Thus p(t) is represented by a constant matrix in the bases {c\{i)}\Li and 
{e K (i)}^ =1 . Let us choose pit) to be the identity in these bases, since this is 
the configuration given by the imbedding. 

The curve git) = (q(t),p(t)) in Isom + (IR 16 ) is given by 

g(t)x = Ax + f(t), 



where 



/ 



Ait) 



cos 2 f 


sin 9 








sin fl 

% fl 


cos fl — 1 
2 „ 











2 


cos^f 








sirr ^ 


sin fl 
2 

















cosf 











sinf 

















1 

















sin fl 
cos fl 2 — 1 


sin 2 f 








cos 2 f 


sinfl 











2 


sin fl 
2 








sin fl 
2 


cos 2 f 

















-sinf 











cosf 





























1 





























cos| 


























sin 5 












— sin 7 
cos 5 



V 



o f 



\ 



0, 



le / 



e e 
\v"'\v 

notes the zero matrix of size mxn and 16 is the identity matrix of size 6x6. 



and r(t) 



-1,0,0,0,0, 



■1,0,0,0,0,0 . Here, mxn de- 
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